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ABSTRACT 

We searched for z = 7.3 Lya emitters (LAEs) behind two gravitationally lensing 
clusters, Abell 2390 and CL 0024, using the Subaru Telescope Suprime-Cam and a 
narrowband filter NB1006 (A^ 1005 nm, FWHM - 21 nm). The combination of the 
fully depleted CCDs of the Suprime-Cam, sensitive to z ~ 7 Lya emission at 1 /xm, 
and the magnification by the lensing clusters can be potentially a powerful tool to 
detect faint distant LAEs. Using the NB1006 and deep optical to mid-infrared images 
of the clusters taken with the Hubble and Spitzer Space Telescopes, we investigated if 
there exist objects consistent with the color of z = 7.3 LAEs behind the clusters. We 
could not detect any LAEs to the unlensed Lya line fiux limit of i^Lya — 6.9 x 10"'^* erg 
s~^ cm^^. Using several z = 7 Lya luminosity functions (LFs) from the literatures, we 
estimated and compared the expected detection numbers of z ^ 7 LAEs in lensing and 
blank field surveys in the case of using an 8m class ground based telescope. Given the 
steep bright-end slope of the LFs, when the detector field-of view (FOV) is comparable 
to the angular extent of a massive lensing cluster, imaging cluster(s) is more efficient 
in detecting z ^ 7 LAEs than imaging a blank field. However, the gain is expected 
to be modest, a factor of two at most and likely much less depending on the adopted 
LFs. The main advantage of lensing-cluster survey, therefore, remains to be the gain 
in depth and not necessarily in detection efficiency. For much larger detectors, the 
lensing effect becomes negligible and the efficiency of LAE detection is proportional 
to the instrumental FOV. We also investigated the NB1006 images of the three z ^ 7 
z-dropout galaxy candidates previously detected in Abell 2390 and found that none 
of them are detected in the NB1006. Two of them are consistent with the predictions 
from the previous studies that they would be at lower redshifts. The other one has a 
photometric redshift of z ~ 7.3, and if we assume that it is at z = 7.3, the unlensed 
Lya line flux would be very faint: FLyQ < 4.4 x 10~^^ erg s~^ cm~^ (Icr upper limit) 
or equivalent width of Wl'^^ < 26A. Its Lya emission might be attenuated by neutral 
hydrogen, as recent studies show that the fraction of Lyman break galaxies displaying 
strong Lya emission is lower at z ~ 7 than at z < 6. 

Key words: cosmology: observations - early universe - galaxies: evolution - galaxies: 
high- redshift. 
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1 INTRODUCTION 

Detecting and characterizing distant galaxies such as Ly- 
man break galaxies (LBGs) and Lya emitters (LAEs) at 
high redshift have been a major method to probe galaxy 
evolution and reionization in the early universe. In addition 
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to the fact that z < 6 LBGs are now routinely observed, 
the number of z > 7 LBG candidates has dramatically in- 
creased recently owning to the great sensitivity and wide 
area view of the Hubble Space Telescope (HST) WFC3 as well 
as new generations of wide near-infrared cameras on ground 
based facilities, and their physical and st atistical properties 
have been intensively investigated ( e.g., 
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Bouwcn s et al. 2010: Bunk er et al.l 2010: Castell ano et al 



2010 : Hickev et al.. 2010: Qesch et al.. .2010: ,Wilkins et al 
201C : iGrazian et al.1 bOlll - lMcLure et all I2OIII : lYan et al 



201l|). On the other hand, th e number of z > 7 LAE samples 
has also iiicreased (lOta et al . 2008, 20T0I: iHiboii et al.ll20ld. 
201ll. l2012l : IShibuva et al.ll2011i : iTilvi et al.ll2010l : lKrug et al.1 



2012r ). and some of LAEs even include sp ectroscopically 
z >, 7 LBG samples llFontana et al 
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Lehnert et al.ll2010':'^Ono et al."2012': 'P eritericci et al 



Schenker et al.l "[2012: VanzcUa ct al. 20 ul). LAEs at 



z > 6 have been used as a probe of reionization, and the 
significant variation in their Lya luminosity function (LF) 
has been observed from 2; ~ 6 out to z ~ 8, which might in- 
dicate either the evolution of LAEs, the possible suppression 
of Lya emission by neutral hydrogen ga s, or th e combination 
of bo t h effects (e.g., Shimasaku et al. 2006; Ka shikawa et al.l 
20061: Iota et al.ll2008l: lOuchi et al.ll2008.: .Hibon et al.l | 2010| : 
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Hibon et al. 2011, ,2012: .Kas hikawa et al. 2011; Krug et al 
2012! : IShibuva et ahlboil ). However, these results still re- 



main tentative for z > 7 LAEs, because the number of can- 
didates detected and studied to date is small. Thus, observ- 
ing more and more z > 7 LAEs and increasing the sample 
size are essential to drawing the statistically more robust 
conclusions about z > 7 Lya LF as well as comparing with 
the z > 7 LBG studies. 

LAEs at z > 7 could be faint and difficult to detect 
due to the possible suppression of Lya emission by neu- 
tral hydrogen, if the universe is not fully reionized. Using 
the magnification by gravitationally lensing clusters is one 
of the effective ways to overcome this problem, as many 
authors have already demonstrated its gre at potential for 
detecting faint high redshift galaxies (e.g., iBouwens et al.l 
20091: iBra dlev et al."2008!: 'Franx et al."l997'; 'Hu et al."2002!; 
Kneibet al. 2004; Richard ct al. 2008; Santos ct al. 2004; 



Stark et al.ll2007l : I Willis et al.l 1200 8). The surface area and 



the flux of an image behind a lensing cluster in the source 
plane are both magnified in the image plane at the time 
of observation. This increases the chance to detect intrin- 
sically very faint objects, if located behind the positions 
of large magnification in the lensing cluster. However, at 
the same time, the larger the magnification factor is, the 
smaller the intrinsic (or "unlensed") surface area is. This 
decreases the chance to detect faint rare objects. In fact, in 
the case of lensing surveys, the unlensed sky area probed 
tends to be small, such as a few to a few tens arcmin^ 
or even smaller, depending on how many clusters are ob- 
served. Usually, the efficiency of galaxy survey is depen- 
dent on the balance between its area and depth as well as 
the slope of galaxy LF. It is suggested that if the slope of 
the LF, — ci(log (/))/ci(log L), is steeper than 1, the gain in 
depth makes up for the loss in are a in the overall number 
of lensed galaxies to be detected ( Broadhurst et al.l 1 19951 : 
iMaizv et al.]|2010l ). iBouwens etHI ^200^ ) actually assessed 



and confirmed this effect by conducting a search for z ^ 7 z- 
dropout galaxies behind 11 massive lensing clusters and find- 
ing 1-4 candidates. In the case of z > 7 LAEs, the bri ght end 
slope of LF tends to be > 1.5. For example, Kashikawa et al.l 
(|201ll) obtained log^* = -3.28 and logL* = 42.76 with 
fixed Q = —1.5 for the Lya LF ba sed on 50 z = 6.6 
LAEs by fitting the ISchechtej l| 19761 ) LF. This best-fit LF 
has — d(log (;/>)/d(log L) > 1.56 at L > L*. Similarly, the 
best-fit Schechter LFs f or 265 and ~ 30 z = 6.6 LAEs from 
lOuchi et all (|2010l ) and IHu et al.l l|2010l) have the slopes of 
— rf(log 0)/d(log Z/) > 1.55 and > 1 .55 at L > L*, respec- 
tively. Moreover. [Hibon et al] (|2010l ) obtained log <j)* = -4.4 
to -3.0 and logL* = 42.5-43.3 with a = -1.5 for the Lya 
LF based on their z — 7.7 LAE candidates, and the slope is 
-d(log0)/d(logL) > 1.55 at L > L*. These LF slopes in- 
dicate that even shallower imaging of lensing clusters could 
be potentially useful for detecting z ~ 7-8 L > L* LAEs. 

Meanwhile, improving detector sensitivity is another 
complementary way to overcome the difficulty in detecting 
z > 7 LAEs. In 2008, the CCDs of the Subaru Telescope 
Suprime-Cam were replaced with the new fully depleted 
CCDs to make it more sensitive to the red wavelength at 
> 0.9 fim and more advantageous for detecting z > 7 galax- 
ies. This had pushed the maximum redshift of galaxies that 
the Suprime-Cam can probe from z = 7 to z = 7.3, and thus 
we made a new narrowband fiher NB1006 (FWHM 21 
nm centered at ~ 1005 nm) designed to target z 7.18- 
7.35 Lya emission. The transmission curve of the NB1006 is 
shown in the Figure[T] Motivated by the potential usefulness 
of the lensing magnification for z > 7 LAE surveys, we con- 
ducted a search for z — 7.3 LAEs behind two lensing clusters 
AbeU 2390 and CL 0024 using the NB1006 combined with 
the fully depleted CCDs on the Suprime-Cam. In this paper, 
we present the results of this search. The paper is organized 
as follows. In Section 2, we describe the imaging observation 
of the lensing clusters and data reduction. Then, the pho- 
tometry and z = 7.3 LAE candidate selection are conducted 
in Section 3. In Section 4, we present the results and discuss 
their implications. We conclude and summarize our results 
in Section 5. Throughout, we adapt a concordance cosmol- 
ogy with {flm,^A, h) = (0.3,0.7,0.7), and AB magnitudes, 
unless otherwise specified. 



2 OBSERVATION AND DATA 

2.1 Optical to Mid-infrared Broadband Data 

The target clusters for our z = 7.3 LAE lensing search are 
AbeU 2390 (z = 0.228) and CL 0024 (z = 0.390). The deep 
optical, near-infrared and mid-infrared broadband images 
taken with the HST and the Spitzer as well as the Subaru 
Telescope are available for AbeU 2390 and CL 0024. These 
data wer e obtained from obse rvations and archives and re- 
duced bv lRichard et al.1 (|2008l ). The optical images of AbeU 
2390 were taken with the HST WFPC2 F555W, F814W and 
ACS F850LP filters (hereafter denoted as V555, Isi4 and zgso, 
respectively), while CL 0024 with the HST ACS F475W, 
F625W, F775W and F850LP fiUers (hereafter 3475, ^625, 
2775 and Z850, respectively). Also, the near and mid-infrared 
images of the both clusters were taken with the HST NIC- 
MOS FllOW, F160W, the Subaru MOIRCS ii'-band and 
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Figure 1. Filter transmission of the Subaru Suprime-Cam nar- 
rowband NB1006 (solid curve) as well as the HST ACS and NIC- 
MOS broadbands zgso, Jiio and Hiao (dashed curves) used for 
our photometry. 

the Spitzer IKAC 3.6 and 4.5 ^im filters (hereafter Jno, 
Hieo, K, [3.6] and [4.5], respectively). The depth of these op - 
tical and infrared images, taken from lRichard et ahl l|2008h . 
are listed in Table 1. 




Figure 2. {Left) The NB1006 image of Abell 2390 cut into the 
3'.2 X 3'.2 size (the NB1006-A2390 image in Section [g^J . (Right) 
The NB1006 image of CL 0024 cut into the 3.'6 X 3.'8 size (the 
NB1006-CL0024 image in Section |221l . The north is up and the 
east to the left. The critical line for a, z = 7.3 source in each 
lensing cluster is shown by the red curve. The areas covered by 
the HST WFPC2 (the cyan line in the NB1006-A2390 image), 
ACS (the cyan line in the NB1006-CL0024 image) and NICMOS 
(the yellow line) are superposed on the images. As described in 
Sections 13.31 and 13.41 we search for z = 7.3 LAEs in the region 
displayed with the cyan line by using the LAE selection criterion 
(1)— (3), and also apply the supplementary LAE selection criterion 
(4) with the NICMOS Jiio band to the objects only in the region 
displayed with the yellow line. 



2.2 Narrowband Observations and Data 
Reduction 

The narrowband NB1006 imaging observations of the lens- 
ing clusters with the Subaru Suprime-Cam were carried out 
at the dark clear nights on 2009 October 15 and 16. At each 
night, they were observed by making use of the interval of 
two hours between the time right after the evening twiUght 
and the time right before the start of our another observa- 
tion, the NB1006 field search for z = 7.3 LA Es. The result 
of the z = 7.3 LAE field search is presented in lShibuva et al.l 
l|201lll paper. We have chosen Abell 2390 and CL 0024 as 
targets because (1) they were observable at that time, (2) 
the deep optical to mid-infrared images taken with the HST 
and th e Spitzer are available for them, and (3) lRichard' et al.l 
(|2008l ) had detected a z-dropout, A2390-z2, whose photo- 
metric redshift has the probability peak at z ~ 7.3, suggest- 
ing that it may have z ~ 7.3 Lya emission. The sky con- 
ditions during observations were photometric with a seeing 
of ~ 0'.'5-0'.'7. By taking the 20 minutes dithered exposure 
frames, we obtained a total of 2.0 and 1.3 hours of imaging 
data for Abell 2390 and CL 0024, respectively. 

We have reduc ed the NB1006 image frames using the 
software SDFRED (|Yagi et al.l l2002l: lOuchi et al.l 12004 ) in 
the same s t andar d manner as in lKashikawa et all (|2004l ) and 
lota et al.1 (|200^ ). The dithered NB1006 exposure frames 
were combined. The seeing sizes of the final combined 
NB1006 images of Abell 2390 and CL 0024 are 0'.'73 and 
0'.'83, respectivel y. The spectrophotometric standard star 
GD71 l|Okelll990l ) was imaged during the observations to cal- 
ibrate the photometric zeropoints of these stacked images, 
which were turned out to be NB1006 — 31.41 and 31.43 mag 
ADU"^ for Abell 2390 and CL 0024, respectively. Abell 2390 
and CL 0024 are located in the center of the large areas of 
NB1006 images taken with the Suprime-Cam (field-of-view 



is 34' X 27'). Hence, we cut the NB1006 images into the 
smaller sizes around Abell 2390 and CL 0024, the similar 
sizes to the effective areas of the WFPC2 and ACS images 
of Abell 2390 and CL 0024, respectively. Hereafter, we call 
them NB1006-A2390 and NB1006-CL0024 images and use 
them for the subsequent photometry and analyses. These 
images are shown in Figure (2] The limiting magnitudes mea- 
sured with a 2'.'0 diameter aperture in the NB1006-A2390 
and NB1006-CL0024 images reached NB1006 = 23.78 and 
23.70 mag at 3a , respectively. The depth of these images 
are also shown in Table 1. 



3 PHOTOMETRY AND CANDIDATE 
SELECTION 

3.1 Photometry 

We conducted photometry to make the NB1006-detected ob- 
ject catalogs. First of all, the astrometry and the pixel scales 
of all the optical to mid-infrared broadband images were 
matched to those of the NB1006-A2390 and NB1006-CL0024 
images. Also, the point spread functions (PSFs) of the 2850 
images were convolved to match those of the NB1006-A2390, 
NB1006-CL0024 and Jno-band images, because we have to 
calculate the zgso — NB1006 and zs5o — Jiio colors by measur- 
ing the 2850, NB1006 and Jno magnitudes using the same 
diameter aperture in order to select z = 7.3 LAE candidates 
(See Section lO)) . 

Then, source detection an d photometry were carrie d out 
with the SExtractor software (|Bertin fc Arnout^ll996l ). The 
pixel size of the NB1006 images taken with the Suprime- 
Cam CCDs is 0'.'202 pixeP^ We considered an area larger 
than five contiguous pixels with a fiux greater than 2a (i.e., 
2x background rms) to be an object. Object detection was 
first made in the NB1006-A2390 and NB1006-CL0024 im- 
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Table 1. Depth of the imaging data 



Abel 2390 




V555 


.f814 


2850 


NB1006 


-'no 


^^160 


K 


[3.6] 


[4.5] 






26.6 


26.2 


26.82 


23.78 


26.27 


26.54 


25.6 


23.9 


23.9 


CL 0024 


9475 


''625 


«775 


2850 


NB1006 


•^110 


^^160 


K 


[3.6] 


[4.5] 




27.81 


27.75 


27.67 


27.28 


23.70 


26.20 


26.60 


25.7 


23.9 


24.1 



NOTES: The depth of NB1006 is 3ct while that of aU the other bands are 5cr taken from iRichard et al.l ijlooi). The following diameter 
apertures were used for the photometry: 0'.'3 for the HST WFPC2 and ACS optical bands, 2'.'0 for the Subaru Suprime-Cam NB1006, 0'.'6 
for the HST NICMOS Jno and Hiao bands, 0'.'5 for the Subaru MOIRCS i^-band and 3'.'0 for the Spitzer IRAC 3.6 and 4.5 bands. 



ages, and then photometry was done in the images of other 
wavebands using the double-imaging mode. The aperture 
magnitudes of detected objects were measured with the 
MAG_APER parameter, and aperture correction s were applied 
to obt ain the total magnitudes. As adopted bv lRichard et al.] 
l|2008l ). we also used the CK'S and (/.'G diameter apertures and 
the corrections of 0.3 and 0.6 mag for the WFPC2/ACS and 
NICMOS images, respectively, for our photometry. These 
aperture corrections were estimated using brig ht isolated un- 
satur ated stars and assuming a point source (|Richard et al.l 
l2008f ). For the measurements of the 2850 — NB1006 and 
.2850 — ^110 colors, we used 2'.'0 and O'.'G diameter apertures. 
We did not apply any aperture correction for NB1006 be- 
cause we did not use total NB1006 magnitude but used only 
the 2" diameter aperture magnitude for measuring the color 
of 2850 — NB1006 and limiting magnitudes of the NB1006- 
A2390 and NB1006-CL0024 images. Finally, by combining 
the photometry in all the wavebands, we constructed the 
object catalogs for Abell 2390 and CL 0024. 




NBIOO6 [AB mag] 

Figure 3. Detection completeness of Abell 2390 and CL 0024 
NB1006 images, calculated for every 0.5 mag bin. The complete- 
ness do not reach 1.0 even for the objects with bright magni- 
tudes. This is because the blended or overlapped objects tend to 
be counted as one object by the SExtractor. 



3.2 Detection Completeness 

What fraction of objects in an image we can reliably detect 
by photometry depends on the magnitudes and overlapping 
of objects. The fraction usually decreases as the magnitudes 
become fainter due to their faintness. Also, the detectabil- 
ity of target objects is affected by their overlapping with 
foreground objects. Especially, this effect could be large in 
the case of lensing clusters where many galaxies are clus- 
tered in a space. To examine what fraction of objects in the 
NB1006-A2390 and NB1006-CL0024 images SExtractor can 
detect or fails to detect to fainter magnitude, we measured 
the detection completeness of our photometry for those im- 
ages. Using the IRAF task starlist, we first created 5,000 
artificial objects with random but uniform spatial and mag- 
nitude distributions, ranging from 20 to 27 mag. Then, we 
spread them over the NB1006-A2390 and NB1006-CL0024 
images, by using the IRAF task mkobject. SExtractor was 
run for source detection in exactly the same way as our 
actual photometry. Finally, we calculated the ratio of the 
number of detected artificial objects to that of created ones 
to obtain the detection completeness. We repeated this pro- 
cedure five times and averaged the obtained completeness. 
The result is shown in Figure |31 The 50% completeness for 
NB1006-A2390 and NB1006-CL0024 images correspond to 
NB1006 ^ 22.5 and ~ 22.75, respectively. 



3.3 LAE Selection Criteria 

Figure [1] shows that the bandpass of NB1006 is located at 
the red side of the bandpass of ACS zgso filter as well as 
in the middle of the bandpass of NICMOS Jno filter. If 
the Lya emission is redshifted in the bandpass of NB1006, 
the LAE is expected to show significant excess in NB1006 
with respect to Z850 and Jno- This characteristic is used 
to isolate z — 7.3 LAEs from other objects. We inves- 
tigated the expected 2350 - NB1006 and Jno - NB1006 
colors for 2 = 7.3 LAEs and derived candidate selection 
criteria. We created model spectra of LAEs by assuming 
the power law continuum fx oc with several different 
slopes /3 — — 3, — 2.5, — 2, — 1.5, — 1, and adding them the 
Lya emission with rest frame equivalent width of Wl'^^ — 
0, 20, 50, 100, 150 and 300A. We did not assume any specific 
line profile or velocity dispersion of Lya emission. Instead, 
we simply added the total line flux value to the spectra at 
1216A. Then, the spectra were redshifted to 2 = 0-8, and 
Lya absorption by intergalactic m edium ( I GM) was applied 
to them, using the prescription of iMadaul (|l995l ). 

Colors of these model LAEs were calculated using their 
redshifted spectra and transmission curves of 2350, NB1006 
and Jno filters and plotted as a function of redshift in Fig- 
ures [J] [S] and [HI For comparison, we also calculated the colors 
of E (el liptical), Sbc, Scd and Im (irregul ar) galaxies by us- 
ing the lColeman. Wu fc WeedmarJ l|l980l ) template spectra 
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Figure 4. zgsg — NB1006 color as a functi on of redshift for various types of galaxies. The colors of E (elliptical), Sbc, Scd and Im 
(irregular) galaxies are calculated by using IColeman. Wu &: WeedmarJ ijlQSOh template spectra. The color of an LAE is calculated 
by using the spectrum fx oc A'' for a variety of UV continuum slopes /3 and Lyo equivalent widths W^'^K Also, colors of L1-L9.5 
and T0.5— T8 type dwarf stars are calculated by using their spectra and plotted (see the text in Section 13.31 for the details of dwarf 
data). (Left) Colors of LBGs (spectra with = 0). Though they do not have Lya emission, LBGs show the narrowband excess of 

2^850 ~ NB1006 ~ 2.3-2.7 at 2 ~ 7-7.3, but amount of excess docs not change much with /3. At these redshifts, the UV continuum is 
detected only at the narrow wavelength range corresponding to the NB1006 passband and the red edge and low sensitivity portion of 
zg50-band. Hence, a difference in /3 makes only a negligible change in the NB1006 and 2g5o fluxes and t hus a mount of excess. (Right) 
Colors of an LBG and LAEs with /3 = —3 are plotted as a representative, as recent study of lOno et al suggests that the 2 ^ 7 

LBGs (2g5o-dropouts) and 2 = 5.7 and z = 6.6 LAEs have the continuum slope of /3 ~ —3 (See Figure|6]for the case of /3 = —2). The 
LAE color clearly shows the narrowband excess at 7.2 < 2 < 7.4. The horizontal and vertical lines indicate the NB1006-excess LAE 
selection criterion we adopt (See Section 13.31 for details). 
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Figure 5. The two color diagram, Jiio — NB1006 versu s 2350 — Jiio, of various types of galaxies. The colors of E (elliptical), Sbc, 
Scd and Im (irregular) galaxies are calculated by using IColeman. Wu &: WeedmanI l|l980l ) template spectra. The color of an LAE 
is calculated by using the spectrum fx oc A'' for a variety of /? and W^Lya' '^^^ symbols indicate the locations of the minimum, 
central and maximum redshifts (2 = 7.18, 7.27 and 7.35 from left to right) of Lya targeted by the NB1006 filter (the redshifts are 
labeled as a example for the /3 = -3, = 150A LAE and the /3 = -3, W£^=^ = 300A LAE in the left and right diagrams, 

respectively). The shaded region shows our LAE selection criteria (3) and (4) (See Section l3.3l l. The solid diagonal line corresponds to 
(2850 — •/no) — (NB1006 — Jiio) = 2850 — NB1006 = 2.3, above which objects satisfy the selection criterion (3). Also, colors of L1-L9.5 
and T0.5— T8 type dwarf stars are calculated by using their spectra and plotted (see the text in Section 13.31 for the details of dwarf 
data). (Left) Colors of LAEs with Wl'^^''^ = 150A are plotted to see the dependency of the colors on /9. Actually, the colors strongly 
depend on the /3. (Right) Colors of an LBG and LAEs with /9 = —3 are plotted as a representative, as recent study of lOno et al.l 1 I2OICII) 
suggests that the 2 ~ 7 LBGs (2850-dropouts) and 2 = 5.7 and 2 = 6.6 LAEs have the continuum slope of ~ —3. The criteria (3) 
and (4) we adopt in Section 3.3 select 2 ~ 7.2-7.4 LAEs with > 50A, z ~ 7.2-7.3 LAEs with VK£™^ = 20-50A and z ~ 7-7.2 

LBGs with VK/<==' = 0. 

Lya 
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Figure 6. The same figures as the right panels in Figures |4] and |5] but for / 3 = —2. We pl o t thes e to see how the colors o f LAEs and 
LBGs change from the case of /9 = —3 because more recent studies such as IWilkins et alj l(20l3) and iDunlop et al. show that 

= —2 for zgso-dropouts. The zgso — Jiio color is more severe and Jiio — NB1006 is more relaxed than the case of /3 = —3. Similarly 
to the case of /3 = —3 (see the right panel of Figure [5]l, the criteria (3) ^850 — NB1006 > 2.3 and (4) zgso — Jiio > 1-5 we adopt in 
Section 3.3 also allow to select z ~ 7.2-7.4 LAEs with VF/°=* > 50A, z ~ 7.2-7.3 LAEs with W/^"* = 20-50A and z ~ 7-7.3 LBGs with 

wrest ^ for 8 = -2. 

Lva " 



and plotted in Figures 21 [S] and [B] Moreover, we also plot the 
colors of L- and T-type dwarf stars, which can be interlopers, 
in Figures |4l [5] and [B] We calculate the colors by using the 
filter transmission curves and the spect r a of L1-L9.5 and 
T0.5-T8 dwarfs ('Burgasser et al.' '2000", '2003*; ' Chiu et al.' 
.2005, 20.06 - Cruz et al. 2003; Fan et al. 2000; G eballe et al 
19961. 12OO1I, (20021; iGolimowslci et al. 1998; Kirkpatricl; et al.l 



1991 1 19991 l200d; iK napp et all 12004; .Leggett et all 



20011 12002a b; M artin et al.l] 19991; iMcLean et al. 
Reid et al. '2000, 2001; Ruiz et al.lll997l;1Schultz et al. 



199S 



2003 



199S 



Strauss et al. 1999; Tsvctanov et al.l 2000l) compiled bv S 



Leggett in her L and T dwarf data archivq 

As clearly seen, an LAE is expected to produce sig- 
nificant flux excess in NB1006 against zgso and Jno at 
around z — 7.3. The model spectra without Lya emission 
(W^Lya = 0), which are considered spectra of LBGs, also 
show some extent of NB1006 excess at 2 ~ 7-7.3. Mean- 
while, elliptical galaxies at z ~ 1-3 also show some excess of 
Z850 - NB1006 ~ 0.9-1.3 due to their 4000A Balmer break. 
Also, L and T dwarfs display excess of zgso — NB1006 ~ 0.5- 
1.9 due to their spectral shapes. Hence, we adopt zgso — 
NB1006 > 2.3 as an NBlOOB-excess criterion that selects 
z ~ 7.2-7.4 LAEs with W^^^ > 50A and z ~ 7.2-7.3 LAEs 
with Wilt = 2O-50A and LBGs (W£°=* = 0) but avoids 
cont aminatio n from the ell ipticals and dwarfs. [Richard et al.l 
(|2008. ) and .Bouwens et al] (j2009. ') conducted surveys of zgso- 
dropout LBGs in Abell 2390 and CL 0024 and detected 
two and one LBG candidates, respectively. Comparison with 
them can distinguish between LAEs and LBGs in our can- 
didate sample. We use the following criteria to select LAEs; 

NB1006 ^ 4a (1) 

Optical bands < 2a and no object is seen visually (2) 
Z850 - NB1006 > 2.3 (3) 

^ http://staff.gemini.edu/~sleggett/LTdata.html. 



Z850 — Jiio > 1.5 (where Jno images are available) (4) 

We set the criterion (1) as our object detection limit. This 
keeps the objects brighter than the ia limiting magnitudes 
(2'.'0 aperture), NB1006 = 23.47 and 23.39 for Abell 2390 
and CL 0024, respectively. The criterion (2) means non- 
detections (< 2a) in V555 and 7gi4 for Abell 2390 and 3475, 
^625 and 1775 for CL 0024 (these are all measured in total 
magnitudes). We also visually inspect the optical bands to 
ensure that no object is seen in each image. This criterion is 
applied since the flux of an LAE shortward of Lya emission 
should be absorbed by IGM and can help reduce the number 
of interlopers such as lower redshift galaxies with other types 
of emission lines (e.g., H/3, [QUI], [Oil] or Ha emitters). 
However, we cannot completely remove all the interlopers 
with this crit erion, as spectroscopy of z = 5.7 and 6.6 LAEs 
conducted bv lShimasaku et al.l (|2006l l and iKashikawa et al.l 
(201 J ) who also applied similar criteria to reduce foreground 
interlopers, for example, identified a few of their LAE can- 
didates as low redshift line emitters (at least 6% and 2% in 
their z — 5.7 and z — 6.6 samples, respectively). 

The criterion (3) is the NB1006 excess against 2:g5o due 
to the z ~ 7.3 Lya emission. The criterion (4) zssq — Jno > 
1.5 is to be applied to the objects at the positions in the 
NB1006 image where the Jno-band image is available to re- 
duce the contamination from low-z line emitters. The colors 
Z850 — NB1006 and zsso — Jiio of the criteria (3) and (4) are 
measured with a 2'.'0 and 0'.'6 apertures. In the calculation 
of the colors, if the zgso and Jno magnitudes are fainter 
than la, they are replaced by the la values. In the right 
panels of Figures [4] and [5] and the both panels in Figure |6l 
we plot colors of LAEs/LBGs with / 3 = —3 and l3 = —2 to 
see the difference as recent study of lOno et al.l (|201G| 1 sug- 
gests that the 2 ~ 7 LBGs (^gso-dropouts) and z = 5.7-6.6 
LAEs have the cont inuum slope of /3 2^ —3 while more re- 
cent s tudies such as IWilkins et al.l (|201lh and iDunlop et al.l 
||2012D show that ^9 = -2 for zgso-dropouts. The zgso — Jiio 
color is more severe and Jno — NB1006 is more relaxed in the 



z = 7.3 LAEs behind Lensing Clusters 7 



case of /3 = —2 than /3 = —3. The criteria (3) and (4) allow 
to select z ~ 7.2-7.4 LAEs with WLy^ > 50A, z ~ 7.2-7.3 
LAEs with Wly^ = 20-50A and z ~ 7-7.2 or 7.3 LBGs with 
Wr^t = for both ^ = -3 and ^ = -2. 



3.4 Candidate Selection 

By applying the LAE selection criteria to the NB1006- 
detected object catalogs, we found 5 and 17 objects sat- 
isfying the criteria in Abell 2390 and CL 0024, respectively. 
Then, we visually inspected their NB1006 images and re- 
moved obviously spurious objects such as tails of saturated 
pixels from bright stars, halos of bright stars and cosmic rays 
with extremely long narrow shapes. As a result, we were left 
with two objects in Abell 2390 and three in CL 0024. The 
two objects in Abell 2390 and two out of the three in CL 0024 
are detected only in NB1006 and not detected in any other 
bands at wavelength longer than NB1006. Two out of them 
are located outside of the Jno and -ffieo images of Abell 2390 
and CL 0024, and we do not know whether they would be 
detected if they were observed in these bands. Hence, there 
is a possibility that these objects might be noises. Mean- 
while, one object in CL 0024 is detected in the K, [3.6] and 
[4.8]. However, this does not necessarily guarantee that the 
detection in NB1006 is real. This object is also located out- 
side of the Jiio and -ffigo images of CL 0024, and we do not 
know whether it would be detected if it were observed in 
these bands. 

Therefore, to check if the five detections in NB1006 are 
real, we visually inspected all the dithered 20 minutes ex- 
posure frames of these objects before they were stacked to 
obtain the final images. We found that at least two out of 
all the 20 minutes frames had cosmic rays with narrow elon- 
gated shapes and high fluxes at the positions of all the five 
objects. These cosmic rays overlapped and were not com- 
pletely removed at the time of the stacking, which caused the 
detections in NB1006. Hence, we conclude that all the five 
objects that satisfied the LAE criteria are not LAEs but the 
residuals of cosmic rays. Usually, cosmic rays can be cleanly 
removed when stacking the dithered exposure frames. How- 
ever, in our case this time, the number of dithered exposure 
frames is not large enough. Hence, if two or more cosmic 
rays accidentally hit the same position in different exposure 
frames and overlap with each other at the time of stacking, 
they cannot be removed completely. Eventually, we could 
not detect any z = 7.3 LAE candidates behind Abell 2390 
and CL 0024. 



3.5 Survey Depth and Volume 

To understand how the depth improved and volume de- 
creased in our survey with the presence of lensing magnifica- 
tion (the trade-off between depth and volume) and to what 
depth and within what area we could not detect z = 7.3 
LAEs, we calculated the unlensed comoving volume sam- 
pled with the NB1006 filter as a function of limiting un- 
lensed Lya luminosity I/Lya in the source plane for each 
cluster. This can be derived by calculating the unlensed 
surface area as a function of magnification jj, (in the unit 
of mag) and converting the area and the NB1006 magni- 
tude unlensed with fi to the comoving volume and the Lya 
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Figure 7. Unlensed comoving volume (surface area) sampled 
with the NB1006 filter brighter than a given unlensed Lya lu- 
minosity (magnitude) in the source plane for each cluster and 
for total of the two clusters. The dotted line shows the total 
volume (area) sampled in the absence of lensing. The limiting 
Lya luminosity shown here is calculated assuming the conversion 
factor of 0.7 between the NB1006 total flux and the Lya flux, 
Fi^ya = 0.7 X Fnb, in the case of z = 7.3 and VK£°=* = 20A (the 
Lya line threshold of our LAE candidate selection criteria). 

luminosity, respectively. The area as a function of /i is cal- 
culated assuming z = 7 . S and using the same lensing model 
used in iRichard et al.l dlOOS'). This calculation is made in 
the field of views (FOVs) of the WFPC2 and ACS for Abell 
2390 and CL 0024, respectively, because the surveyed areas 
in the image planes of Abell 2390 and CL 0024 are limited 
by the WFPC2 and ACS images used for the LAE selec- 
tion criterion Q (null detections in the bands blueward of 
Lya) combined with the criterion (3). The conversion of the 
area to the volume is done by considering the comoving dis- 
tance along the line of sight corresponding to the bandwidth 
(FWHM) of the NB1006 filter. 

When calculating the Lya luminosity from the NB1006 
magnitude, we used the conversion factor of 0.7. That is, we 
consider that 70% of the flux corresponding to the NB1006 
magnitude comes from Lya emission. This conversion fac- 
tor is derived as follows. We assume that the NB1006 flux 
-Fnb consists of contributions from Lya and UV contirmum 
fluxes, Fi^ya and Fcont- 

-Fnb = Fi^ya + -Fcont (5) 

The Lya flux is related to the UV continuum flux density 
/A,cont with observed Lya equivalent width W^yJ. 

Fhya = VKLyQ,fA,cont (6) 

On the other hand, the /A,cont is the UV continuum flux 
divided by the wavelength from Lya emission to the edge of 
the NB1006 fllter bandpass (Anb" - A^!;;";), 

/A,cont = -fcont/(ANB'' ~ Aj^ya ) (7) 

where A}?!'' = IOI6OA, and we assume that Lya emission is 
at 2 = 7.3, i.e., A^!;;^^ = (1 + z)1216A = 10092.8A. Solving 
these equations ©-([T]) for Fi^ya, we have 

/ \ max \ obs \ 

F^... = w(l + ^^^%g^^) (8) 
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Figure 8. Distribution of magnification fac tors ii (in magn i tudes ) 
calculated by using the lensing models of iRichard et all l|2008l ) 
and assuming a point source at z = 7.3. The horizontal axis is 
/i in 0.1 mag bin while the vertical axis shows unlensed surface 



Because our LAE selection criterion (3) corresponds to the 
rest frame Lya equivalent width threshold of W][y2 ^ 20 A, 
we have W£^l = WUti^ + z) ^ 166 A and accordingly 
Fhya ^ 0.7Fnb. The ^Lya ~ 0.7Fnb or equivalently I/Lyo, ~ 
0.7I/NB in luminosity corresponds to the limiting luminosity 
in our lensing survey and the volume should be calculated 
down to this limit. Figure [7] shows the volume sampled with 
the NB1006 as a function of limiting Lya luminosity in the 
source plane for each cluster and the sum of the two clusters. 
The trade-off of increase in depth and reduction in volume 
is clearly seen. In Figure[Sl we show the distribution of mag- 
nification factors fj, (unlensed area versus /j, in 0.1 mag bin) 
in each cluster. The /i ranges from to 2.5 mag and the 
total unlensed area (sum of Abell 2390 and CL 0024) in the 
source plane is 8.4 arcmin'^, a factor of 2 reduction from 
the lensed area in the image plane (i.e., sum of the FOVs 
of the WFPC2 and ACS; see Figure 01. The unlensed 4a 
limiting magnitudes in the fi = 1.3 (median magnification) 
and fj, = 2.5 areas are NB1006 = 24.7 and 25.9. These corre- 
spond to the Lya line flux limits of i^Lya — 2.1 x lO"'^'^ and 
6.9 X 10"^® erg s"^ cm'^ or Lljc ^ 1.3 x 10^^ and 4.4 x 10*^ 
erg s~^, respectively. Therefore, to the limits varying from 
LLya ^ 4.4 x 10"*^ (/X = 0) to 4.4 x lO'^^erg s"^ (as shown in 
Figure [7]) within the total area of 8.4 arcmin'^, we could not 
detect any z = 7.3 LAEs. 



4 DISCUSSION 

4.1 Comparison with Blank Field Surveys 

To understand the meaning of the null detection of 2 = 
7.3 LAEs in our lensing survey, we estimate the expected 
detection numbers of z ~ 7 LAEs in lensing surveys and 
compare them with those in blank field surveys in the case 
of spending the same amount of observing time with an 8m 
class ground based telescope, a camera with CCDs sensitive 
to z ~ 7 Lya (red-sensitive CCDs), and a narrowband filter. 
With this estimation and comparison, we can see if the null 
detection in our lensing survey was an anticipated result 



and know what is the most efficient way to detect 2 ~ 7 
LAEs. We consider a model case that an 8m class telescope 
= Subaru Telescope, red-sensitive CCDs = Suprime-Cam's 
fully depleted CCDs, a narrowband filter = NB973 (Ac — 
9755A, FWHM = 200A) and a lensing cluster = CL 0024. 
The NB973 i s the filter used bv lOta et all ||2008| . |2010D and 
iHibon et all (|2012l ) for their 2 = 7 LAE surveys. 

We first calculate the expected detection number of 
LAEs in lensing and field surveys by using Lya LFs derived 
from the several differ ent narrowband z = 7 LAE blank field 
sur veys conducted bv|Qta et al.1 (|2010l '). iHibon et all (1201 ll ) 
and iHibon et al.1 (|2012l ) (hereafter OlO, Hll and H12) and 
then consider the case for the narrowband NB1006 z = 7.3 
LAE lensing survey (the current study). We use the three 
different 2 = 7 LFs for the calculation since we do not have 
any z — 7.3 LFs and we also want to see how the expected 
number of LAEs changes with different LFs. 

We compare four cases (see also Table [2|: Case (1) a 
survey whose one pointing covers the area comparable to 
the HST ACS field of view (FOV ~ 11.3 arcmin^ and corre- 
sponds to the area just covering the size of a lensing cluster 
CL 0024) and observing many clusters (1, 2, 5, 10, 50 and 
100 X CL 0024), Case (2) a lensing survey covering an area 
with the HST ACS FOV and observing only one cluster (CL 
0024) to several different depths. Case (3) a survey observ- 
ing one blank field covering an area with the ACS FOV 
to several different depths and Case (4) a survey observing 
one blank field covering an area comparable to the Subaru 
Suprime-Cam FOV to several different depths. For the Case 
(1), we consider integration time per cluster of 2 hours as 
a model case, because spending short integration time for 
each cluster by counting on the magnification and observing 
many clusters are considered a practical observation strat- 
egy with a ground based telescope. For the Cases (2)-(4), 
we consider integration times of 2, 4, 10, 20, 100 and 200 
hours, which corresponds to 2 hours x 1, 2, 5, 10, 50 and 
100 clusters, to compare with the result of the Case (1). 



4.1.1 The Expected Detection Number of z • 
Lensing Surveys 



7 LAEs m 



We make the calculation of this with the following 5 steps. 
(1) First, we convert the magnification /i distribution of CL 
0024 in Figure|8]to unlensed comoving volume V{L) sampled 
down to each unlensed limiting Lya luminosity (LLya) bin. 
In this procedure, we convert each /j, to unlensed 5a NB973 
limiting magnitude with m\f^ + fi and then to the corre- 
sponding I/Lya , where the m\Z" is lensed 5a NB973 limiting 
magnitudes reached by the integration times we consider. 
We use the conversion factor used by OlO to convert NB973 
magnitudes to LLyo, for the OlO case and the conversion 
factor used by Hll for the Hll and H12 cases. (2) Next, we 
calculate cumulative LAE number density n(>L) down to 
each LLya bin by linearly extrapolating the 2 = 7 Lya LF 
(log n(>L) vs. logLLya) obtained by OlO. This LF has a 
slope of — d(log n)/d(log LLya) — 1.5. Meanwhile, to obtain 
n(>L )'s for Hll and H12, we integrate the best-fit [Schechteil 
l|l97(j ) 2 = 7 Lya LFs derived by them to each LLya bin (See 
the footnote of Table [2] for more details). (3) By multiplying 
each V{L) obtained in (1) by n{>L) obtained in (2), we cal- 
culate the number of LAEs N{L) = V{L)n{>L) expected 
to be detected down to each I/Lya bin. (4) Summing up all 
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Table 2. Comparison of the expected detection number of 2 ~ 7 LAEs in lensing and blank field surveys using an 8m class ground-based 
telescope, red-sensitive CCDs and a narrowband filter. 



Case (1) Case (2) Case (3) Case (4) 



Obs. Observing many lensing clusters Observing 1 lensing cluster Observing 1 blank field Observing 1 blank field 



Time" 


Number of 


#LAE; 


3 in ACS FOV 


#LAEs in ACS FOV 


#LAEs 


in ACS FOV= 


#LAEs 


in S-Cam FOV^ 


(hour) 


Clusters* 


010 


Hll 


H12 


010 


Hll 


H12 


010 


Hll 


H12 


010 


Hll 


H12 


2 


1 


0.02 


0.22 


0.03 


0.02 


0.22 


0.03 


0.01 


0.02 


0.01 


0.86 


1.5 


0.47 


4 


2 


0.04 


0.45 


0.06 


0.04 


0.61 


0.07 


0.02 


0.14 


0.03 


1.73 


9.9 


2.0 


10 


5 


0.11 


1.1 


0.15 


0.11 


2.1 


0.21 


0.06 


0.87 


0.12 


4.38 


63.7 


8.7 


20 


10 


0.21 


2.2 


0.30 


0.21 


4.4 


0.40 


0.12 


2.5 


0.28 


8.83 


183 


20.3 


100 


50 


1.1 


11.2 


1.5 


1.1 


17.4 


1.3 


0.62 


14.4 


1.2 


45.0 


1048 


85.9 


200 


100 


2.1 


22.4 


3.0 


2.1 


27.3 


1.9 


1.2 


24.9 


1.9 


90.7 


1819 


138 


#LAEe 


i per hour'' 


0.01 


0.11 


0.02 


0.01^ 


0.22-f 


0.02^ 


0.01^ 


0.13-^ 


0.01^ 


0.44^ 


9.17-'' 


1.02^ 


Hours 


per LAE^ 


95 


8.9 


66 


96-f 


4.6-f 


50-f 


165^ 


8.0^ 


72f 


2f 


0.1-'' 


1.0-f 



NOTE: This is a model case that an 8m class telescope = Subaru, red-sensitive CCDs = Suprime-Cam fully depleted CCDs and a 
narrowband filter = NB973 (Ac = 9755A, FWHM = 200A). "Total observing time in the unit of hour. In the Case (1), we considered 
the integration time per cluster of 2 hours, and hence total observing time = 2 hours X the number of clusters. In the Cases (2)-(4), the 
total observing time means integration time of one cluster or one blank field. 'The number of lensing clusters observed for the Case (1). 
'^The number of LAEs expected to be detected in the field of view (FOV) of th e ACS and the Suprime-Cam estimated by extrap olating 
the z = 7 Lyo LF derived by OlO (Ota ct al. 2010.) and integrating the best-fit ISchechteJ ([w76) Lya LFs derived by Hll ( Hibo n et all 
l201lh and H12 (Hi bon et al.ll2012f l. In the case of H12, we use their "common sample LF" listed in their Table 4. Because logL* = 43.71 
listed in their Table 4 was found to be a typo, we use the correct one, logL* = 42.71 (P. Hibon 2012, private communication) for our 
calculation. ''The number of LAEs expected to be detected in one hour. ''The time (in hour) necessary to detect one LAE. ■''The values 
in the case of integration time of 20 hours (a practical observing time for an 8m class ground-based telescope) are listed because the 
expected detection number of LAEs varies as slope of the Lya LF changes from the bright end to the survey depths, and the survey 
depth depends on the integration time. 



the N{L) in all the I/Lya bin, we obtain the number of LAEs 
expected to be detected behind one lensing cluster (in this 
case, CL 0024). (5) The number of LAEs expected to be 
detected in surveying more than one clusters is calculated 
by simply multiplying the number of LAEs obtained in (4) 
by the number of clusters in consideration. 

4.1.2 The Expected Detection Number 0/2: ~ 7 LAEs in 
Blank Field Surveys 

We make the calculation of this with the following 3 steps. 
(1) Linearly extrapolating the OlO's z = 7 Lya LF or 
integrating the Hll's and H12's best-fit Schechter z = 7 
Lya LFs, we calculate the cumulative LAE number densi- 
ties n(>L)'s down to the limiting LLya's corresponding to 
5a limiting NB973 magnitudes reached by the integration 
times we consider. (2) Multiplying these n(>I/)'s by the vol- 
ume sampled by the Suprime-Cam -I- NB973 filter, 3 x 10^ 
Mpc'^ (corresponding to an area of 824 arcmin^) taken from 
010, we obtain the number of LAEs expected to be detected 
in the blank field survey with the Suprime-Cam FOV. (3) 
Then, we calculate the expected detection number of LAEs 
in the blank field survey with the same sampled area as the 
lensing surveys. In the case of CL 0024, the area corresponds 
to the FOV of ACS ~ 11.3 arcmin^. Hence, we obtain the 
number of LAEs by multiplying the number of LAEs in the 
Suprime-Cam blank field survey by a factor of 11.3/824. 

4.1.3 Comparison of Lensing and Field LAE Surveys 

The result of above calculations in Sections 14.1.11 and 14.1.21 
are presented in Table[2l First, we compare the Cases (l)-(3) 
for each LF. In the case of the OlO LF, the survey efficiencies 



of the Cases (1) and (2) are the same and better than that 
of the Case (3). This is expected because the OlO LF has 
the bright end slope of — d(log n)/d(log LLya) — L5 > 1 and 
we simply linearly extrapolate this LF to the fainter lim- 
its corresponding to the integration times to calculate the 
number of LAEs. In the case of the Hll LF, the Case (2) is 
more efficient than the Cases (1) and (3), and the Case (1) is 
more efficient than the Case (3) until the Case (1) observes 
10 clusters and the Case (3) integrates for 20 hours; there- 
after, the trend is reversed for the longer observing times. 
Meanwhile, in the case of the H12 LF, the Case (2) is more 
efficient than the Cases (1) and (3) until the observing time 
reaches 20 hours; thereafter, the efficiencies are comparable 
between the Cases (2) and (3), and that of the Case (1) is 
better than them. 

In summary, for the observing times shorter than 20 
hours, the Case (2) where we concentrate on one lensing 
cluster is more efficient than observing many clusters with 
short integration (Case (1)) or observing one blank field with 
an area comparable to size of the lensing cluster (Case (3)), 
but the trend for the longer observing times varies with LFs. 
All the three LFs considered here have their observed data 
points only at the bright end (i.e., a few z = 7 LAE candi- 
dates with Li^ya ~ 0.9-1 X 10^^ erg s~^), and we extrapolate 
such LFs linearly or with the Schechter function. This re- 
sults in larger uncertainties in the slopes of the LFs at the 
fainter luminosity limits reached by observing times we con- 
sider. Hence, the results in Table[5]for the shorter observing 
times refiect the observed data points of each LF better than 
those for longer observing times. In reality, spending more 
than 20 hours for an observation with a ground based 8m 
telescope is very expensive and not practical. Hence, con- 
centrating on one lensing cluster seem to be most efficient 
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for detecting z = 1 LAEs, but two (OlO and H12) out of 
the three LFs imply that even spending 20 hours, we de- 
tect 0.2-0.4 LAEs (null detection) while the other one LP 
(Hll) predicts 4.4 LAEs. 

Next, we compare this result with the Case (4) where 
we observe one blank field whose area is much larger than 
the size of a lensing cluster (i.e., the Suprime-Cam FOV in 
the current case) . In the Case (4) , the three LFs predict the 
detections of 8.8-183 LAEs for 20 hours of observing time, 
much larger than 0.2-4.4 LAEs in the Case (2) considered 
above. Also, for any observing times and all the three LFs, 
the Case (4) is much more efficient than the Cases (l)-(3) 
in detecting z = 7 LAEs. Hence, in conclusion, observing 
sufficiently large sky area compared to the size of lensing 
clusters is the most efficient way to detect z ^ 7 LAEs with 
an 8m class ground based telescope, a red-sensitive CCDs 
and a narrowband filter of 200A width. 

However, three things should be noted. First, we used 
the magnification (fi) distribution of CL 0024 as a model 
case for the calculation of the expected number of LAEs. 
The jj, distribution varies from cluster to cluster, and thus 
the result presented in Tableland the text above might not 
necessarily represent the most common or the most proba- 
ble case for lensing surveys of « ~ 7 LAEs. Second, as al- 
ready mentioned, the uncertainty is large for the results with 
longer observing times due to the large uncertainty in the 
faint end slopes of the LFs. Finally, all the LFs considered 
here consist of photometric LAE candidates. If they include 
some contaminations, the real slopes of the LFs are different 
and thus so does the result in Table [2] The difference in the 
expected number of LAEs predicted by the three LFs could 
come from the difference in the degree of contaminations. 



4-1-4 Expected Detection Number of z = 7.3 LAEs m the 
Present Lenstng Survey 

We also make the similar calculation for the current z = 
7.3 LAE lensing survey (2.0 hours for A2390 -I- 1.3 hours 
for CL 0024; lensed depth of NB1006 ~ 23.4 at 4cr) with 
the narrowband NB1006 to obtain the expected detection 
number of LAEs by following the same procedures in Section 
14.1.11 and assuming that the z = 7.3 Lya LF is similar to 
z = 7 Lya LF. We find that the expected detection number 
of LAEs is 0.0062, 0.081 and 0.0096 if we use the OlO, Hll 
and H12 Lya LFs. This is consistent with the null detection 
of LAEs in our survey. Since the Suprime-Cam sensitivity 
at the wavelength of 2 = 7.3 Lya is lower than that at z = 7 
Lya, detecting the z = 7.3 LAEs in a lensing survey is more 
difficult. 



4.2 Comparison with z-dropouts Previously 
Detected in Abell 2390 

Two searches for z >^ 7 zgso-dropout galaxies have been 
conducted in t he len sing clusters Abell 2390 and CL 0024. 
iRichard et"all (j2008h detected two zs s o-drop outs, A2390-zl 
and A2390-z2, while iBouwens et al.l ()2009l ) detected one, 
A2390-zDl. Interestingly, one out of the three, the A2390- 
z2, has a photometric redshift whose prob ability peaks at 
z ~ 7.3 (see Figure 7 of lRichard et aLlbOOSj ). Also, the mag- 
nification by the lensing cluster at the position of this object 



is ^ = 1.8 mag (see Table 3 of lRichard et al.ll2008l ). Hence, if 
this object is really nt z ^ 7.3 and has Lya emission, its line 
fiux could be boosted by the magnification and detected in 
NB1006. However, we found no object detected in NB1006 
at the position of the A2390-z2, as seen in Figure [9] There 
are three possible interpretations. (1) The A2390-z2 is not at 
z ~ 7.3. (2) It is at 2; ~ 7.3 but does not show Lya emission. 
(3) It is at z 7.3 and have Lya emission, but the emission 
is fainter than the NB1006 image detection limit. We cannot 
tell which is most likely from the current data alone. If the 
A2390-z2 is at z ~ 7.3, it might be possible that its Lya 
emission is suppressed by neutral hydrogen^ as recent stud- 



ies llStark et al 2010. 20 yl; lOno et al.ir2012l : (Pentericci et all 
I2OIII : Schenker et al, 2013 ) suggest that the fraction of Lya 



emitting LBGs would be smaller at 2: > 6-7 due to such 
an effect. If we assume that the A2390-z2 is at z ~ 7.3, we 
estimate the unlensed upper limit on the Lya line fiux to be 
Fijya < 4.4 X lO"'^* erg s~^ cm~^, by directly converting the 
Icr limiting magnitude of NB1006 = 24.98 plus /x = 1.8 mag 
into the fiux. As seen in Figure [9] the A2390 -z2 is detected 
in H\Qo with iifieo = 25.8 mag (see Table 3 of lRichard et al.l 
[2OOS'). This magnitude can be converted to the continuum 
fiux density of /A,cont ~ 2.0 x 10"^° erg s"^ cm"^ A. Using 
this and equation (jSJ and assuming z = 7.3, the upper limit 
on Fi,ya corresponds to W^y^ < 26A. 

We also checked the NB1006 images of the other two 
objects, A2390-zl and A2390-zDl, and confirmed no detec- 
tions for the both cases (see Figure |9]). The A2390-zl has a 
photometric redshift whose prob ability peaks at z ~ 6.4 
(Figure 7 of iRichard et al.l |2008| ). and the non-detection 
in NB1006 is consis tent with this. On the other hand, 
iBouwens et~aLl (|2009l ) selected the A2390-zDl as a "weak" 
zssQ-dropout based on the same selection criteria adopted by 
IBouwens et al.l (|2008l ) . Comparison of the colors and photo- 
metric errors of the A2390-zDl, zgso — Juo = 1-1 ± 0.8 an d 
Jiio - H160 = 0.8 ± 0.2 (Table 2 of IBouwens et al.ll2009l ). 
with the zgso-dropout selectio n criteria in the zgso — ^iio vs. 
Jiio — Him plane (Figure 2 of IBouwens et al.ll2008l ) suggests 
that this object is at the border of the 2850-dropout crite- 
ria and can be either at 2: ~ 6.5-7.0 (inside the criteria) or 
z < 6.5 (outside the criteria), depending on the combina- 
tion of the photometric errors in the two colors. Since these 
redshifts are lower than 2 = 7.18, the minimum redshift of 
Lya the NB1006 can detect, non-detection in NB1006 is a 
consistent result. 



5 CONCLUSION 

We imaged the two lensing clusters, Abell 2390 and CL 
0024, using the narrowband NB1006 (central wavelength 
~ 1005 nm and FWHM ~ 21 nm) with the Subaru Tele- 
scope Suprime-Cam to search for gravitationally lensed 
2 = 7.3 LAEs. The fully depleted CCDs of the Suprime- 
Cam have the sensitivity to 2 ~ 7 Lya emission at ~ 1 
^m. Also, recent studies of 2 ~ 6.6-7.7 LAEs show that 
the Lya LFs at 2 ~ 7 have the steep bright end slope of 
— d(log0)/ci(logL) ~ 1.5 > 1 at 1/ > L*, suggesting that 
even shallower imaging of lensing clusters could be poten- 
tially useful for detecting 2 > 7 bright LAEs with the gain 
in depth by lensing magnification possibly making up for 
the loss in area in the overall number of lensed LAEs to be 
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Figure 9. The 4" X 4" broadband and NB1006 i mages of ^ssp- 
dropo uts previously detected in Abell 2390 by [Richa rd et alj 
1I2OO8I ') (A2390-Z1 and A2390-z2) and iBouwens et alJ ll2009f) 
(A2390-zDl). The A2390-z2 has a photometric redshift whose 
probability peaks at ^ ~ 7.3. All the three zgso-dropouts are not 
detected in NB1006. 



detected. Hence, combination of the Suprime-Cam CCDs 
and the lensing magnification was considered to be a pow- 
erful strategy to detect z > 7 LAEs. The magnification in 
Abell 2390 and CL 0024 varies from ^ = to 2.5 mag. We 
reached the lensed depth of NB1006 ~ 23.4 at 4a in both 
lensing cluster images. Where fj, — 1.3 (median) and /i = 2.5, 
this corresponds to the unlensed depths of NB1006 ~ 24.7 
and 25.9 or Lya line fiux limits of -pLya — 2.1 x 10^ and 
6.9 X 10~^* erg s~^ cm"'^, respectively. 

Performing the photometry in the NB1006 images as 
well as deep optical to mid-infrared images of the clusters 
taken with the HST, Subaru and Spitzer, we constructed 
the NB1006-detected object catalog and applied the selec- 
tion criteria of z = 7.3 LAEs. We could not detect any LAEs 
to our survey limit. Using the z = 7 Lya LFs determined 
by three different previous surveys (OlO, Hll and H12), we 
estimated and compared the expected detection numbers of 
z ^ 7 LAEs in lensing and blank field surveys with an 8m 
class ground based telescope, red-sensitive CCDs and a nar- 
rowband filter. We found that a blank field survey covering 
an area sufficiently larger than lensing clusters is more ef- 
ficient in finding z ~ 7 LAEs than both a lensing survey 
observing many clusters with shallow imaging and a lensing 
survey imaging one cluster to a deeper luminosity limit, ex- 
pected from the bright end slopes of the three z ~ 7 Lya 
LFs. 

We also investigated the NB1006 images of the three 
z ^ 7 z-dropouts previously detected in Abell 2390 and 
found that none of them are detected in the NB1006. Two 
of them are consistent with the predictions from the previous 
studies that they would be at lower redshifts. The other one 
has a photometric redshift of 2 ~ 7.3, and if we assume that 
it is at 2 = 7.3, the unlensed Lya line flux would be very 
faint: F^ya < 4.4 x 10"^* erg s~^ cm~^ {la upper limit) 
or correspondingly Wly^^ < 26 A. Its Lya emission might be 
possibly suppressed by neutral hydrogen, as recent studies 
suggest the smaller fraction of Lya emission among LBGs 
at 2 > 7 than those at lower redshifts due to such an effect. 
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